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ABSTRACT 

The Rapidly Deployable Stable Platform (RDSP) is a novel 
vessel designed to be a reconfigurable, stable at-sea platform. 
It consists of a detachable catamaran and spar, performing mis- 
sions with the spar extending vertically below the catamaran and 
hoisting it completely out of the water. Multiple thrusters located 
along the spar allow it to be actively controlled in this configura- 
tion. A controller is presented in this work that uses an adaptive 
feedback algorithm in conjunction with Direct Adaptive Distur- 
bance Rejection (DADR) to mitigate persistent, vortex-induced 
disturbances. Given the frequency of a disturbance, the nominal 
DADR scheme adaptively compensates for its unknown ampli- 
tude and phase. This algorithm is extended to adapt to a dis- 
turbance frequency that is only coarsely known by including a 
Phase Locked Loop (PLL). The PLL improves the frequency es- 
timate on-line, allowing the modified controller to reduce vortex- 
induced motions by more than 95% using achievable thrust in- 
puts. 


LIST OF ACRONYMS 

DADR 

Direct Adaptive Disturbance Rejection 
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STD 
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Figure 1 . THE RDSP’S TWO OPERATIONAL MODES. 


INTRODUCTION 

Quasi-steady currents around a spar can result in 
disturbance-generating vortices. These are shed at relatively low 
frequencies [1], and the resulting forces are nearly periodic and 
somewhat theoretically predictable [2], These frequencies can 
be estimated empirically, but cannot be precisely determined us- 
ing existing analytic means. Their effect has been studied for 
offshore structures in [3], and specifically for cylinders or spar 
platforms in [2,4, 5]. If the frequency of vortex shedding is 
in resonance with the natural frequency of the moored struc- 
ture that produces it, large amplitude Vortex-Induced Vibrations 
(VIVs) can develop, introducing large stresses on the member 
and mooring lines [1]. The use of active control for moored spar 
platforms to mitigate VIVs was introduced in [2]. For a spar 
with a length of 198 m and a diameter of 38 m, sixteen 3000 hp 
thrusters were recommended at a total cost of between 7 and 8 
million dollars [2]. This proposed active control system resulted 
in a significant reduction in installation cost when compared to 
the typical passive VIV-control devices (such as helical strakes). 

This paper will consider vortex-induced disturbances on an 
un-moored powered spar platform called a Rapidly Deployable 
Stable Platform (RDSP). The RDSP is currently under develop- 
ment to create a mobile yet stable offshore platform for multiple 
applications. It consists of two elements: i) a catamaran and if) a 
long, detachable spar buoy, as depicted in Fig. 1. The RDSP can 
be configured for two modes of operation: Horizontal Transit 
Mode (HTM) and Vertical Operations Mode (VOM). The ma- 
jority of the potential missions for the RDSP will be conducted 
when it is operating in VOM, where spar extends vertically below 
the catamaran and lifts it completely out of the water (Fig. 1). 
This spar configuration presents a small cross-sectional area to 
passing waves, thereby improving its seakeeping properties. The 


spar component of the RDSP is equipped with 12 thrusters for 
control in the VOM. These thrusters allow the RDSP to transit at 
slow speeds, dynamically position, and mitigate disturbances. A 
one tenth scale model of a RDSP (1 1 .4 m spar length) has been 
designed and constructed at Florida Atlantic University (FAU), 
complete with a sensor suite that is capable of providing all of 
the states of the RDSP [6]. Preliminary findings for a full scale 
RDSP with a 10 m diameter reveal vortex-induced oscillations 
with displacement amplitudes of approximately 1 m [7]. This is 
significant when dynamically positioning near other vessels for 
cargo transfer applications. Mitigating these disturbances for un- 
moored applications is the focus of this paper. 

An adaptive control approach that uses Direct Adaptive Dis- 
turbance Rejection (DADR) to mitigate vortex-induced distur- 
bances will be taken in this paper. The disturbance rejection 
algorithm will be appended with a Phase Locked Loop (PLL) 
for on-line disturbance frequency estimation. This is suitable for 
mitigating vortex-induced RDSP motions since they are associ- 
ated with eddy shedding frequencies that can be estimated but 
cannot be precisely determined. The controller will take advan- 
tage of the known and estimated disturbance properties (wave- 
form and initial frequency guess) and adaptively compensate for 
the unknown or estimated components (amplitude, phase and fre- 
quency). Results will be simulated using a 6 Degree Of Freedom 
(DOF) computer simulation of this powered ocean platform de- 
veloped in [8,9]. 

The development of the adaptive control and disturbance re- 
jection algorithm used here may be found in [10], prior to aug- 
menting it with a PLL. This was previously applied to the RDSP 
without the inclusion of frequency synthesis for wave and load 
disturbance rejection in [1 1-13]. The concept of augmenting the 
algorithm with a PLL in order to reject disturbances at coarsely 
known frequencies was considered in [7, 14], and applied to a 
simplified 1 -DOF RDSP simulation. This paper extends the work 
of [7] by providing a more detailed description and analysis of 
the PLL, and applying it to a more realistic and complex simu- 
lation of the RDSP that utilizes all 6-DOF (rather than 1-DOF). 
Furthermore, multiple sensor and actuator configurations will be 
considered. 

The paper will be organized in the following manner. First, 
we will begin with a summary of the 6 DOF model for the RDSP, 
including the current-driven vortex forces. Next, the adaptive 
control and disturbance rejection framework will be introduced, 
including a modification for disturbances at frequencies that are 
not precisely known. Simulation results are shown for two RDSP 
scenarios: (i) a single output at the top of the spar controlled us- 
ing a single force located at the center thruster (ii) control actu- 
ation via force and moment generated at the center thruster with 
two collocated outputs. Finally, some concluding remarks will 
be made. 
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NUMERICAL MODEL OF THE 1/10™ SCALE RDSP 
PROTOTYPE 

The control systems presented in this paper are developed 
and tuned to control the \jW h scale RDSP prototype designed 
and constructed at Florida Atlantic University. The initial design 
of this system is presented by [1 5] and the upgraded prototype is 
presented in [8] and briefly reviewed here. This prototype has a 
cylindrical spar section with a length of 1 1 .4 m and a diameter of 
0.84 m (Fig. 1). The spar is attached via a hinge to a catamaran 
with pontoons that both have lengths of 4. 1 m and diameters of 
0.42 m. The catamaran has a dry mass of 361 kg while the spar 
has a total dry mass of 1 , 437 kg, can hold 1 , 3 1 7 kg of ballast wa- 
ter, and has 4, 79 1 kg of entrained sea water, not including ballast 
water, when fully submerged. When operating in VOM (Fig. 1 
(b)) the catamaran and spar are rigidly locked together forming 
a single nearly rigid body that has a total wet mass including en- 
trained sea water of up to 7,906 kg. When configured in HTM 
(Fig. 1 (a)), the ballast tanks are only partially filled with water 
and the entrained sea water is reduced, thereby reducing the total 
wet mass of the RDSP. 

The 1/1 0th scale RDSP prototype has 12 thrusters that are 
located on the spar to allow the RDSP to actively control its 
position and orientation (Fig. 1). These thrusters are automat- 
ically controlled using a computer system that links state feed- 
back from a comprehensive sensor suite including an IMU, dual 
antenna DGPS compass, and tilt sensor to thruster inputs. Each 
of these thrusters can individually produce a maximum forward 
thrust of 320 N and a maximum reverse thrust of 170 N. This 
is much greater than the maximum scaled thrust of 47 N per 
thruster estimated from the proposed generator system suggested 
for the full scale system [8,9]. Not only is this system control- 
lable using thrusters, but it also has an active ballasting system 
that can be used to change the spar’s inertial properties, raise and 
lower the catamaran, and transfer the RDSP prototype between 
HTM and VOM. 

6-DOF Mathematical Model 

A 6-DOF mathematical model and computer simulation of 
the 1/10^ scale RDSP prototype has been developed. This sim- 
ulation models drag, added mass, wind, vortex, and hydrostatic 
forces on the RDSP. 

Kinematics When the catamaran and spar are rigidly 
connected in VOM (Fig. 1 (b)), the RDSP can be modeled as 
a rigid body with its states defined using two coordinate systems. 
These coordinate systems are the earth fixed coordinate system, 
Je, denoted by capital variables and body fixed coordinate sys- 
tem, denoted by lowercase variables. The earth fixed coordi- 
nate system, /Fe, has its origin at the surface of the mean seawater 
with the X-axis pointing North, 7-axis pointing East, and Z-axis 
pointing down. The body fixed coordinate system, J^ 9 is fixed to 


the RDSP with its origin where the spar and the catamaran attach 
and the x-axis pointing towards the bow of the catamaran, they- 
axis pointing towards starboard, and z-axis pointing downwards 
along the spar; in accordance with the right hand rule (Fig. 1 (b)). 
An additional coordinate system is used by the MIMO controller, 
Jc 9 has its origin midway between the center four thrusters with 
its axes aligned with those of Je and defined by xc, yc > and zq 
(Fig. 1 (b)). 


Drag and Added Mass Hydrodynamic drag and added 
mass forces are modeled using Morison’s equation, which ac- 
counts for the interactions between wave, current, and body 
(catamaran and spar) motions [16], 

/= }^C&pAu\u\ +C a puV, (1) 


where p is the density of sea water, ( 7 ) indicates that the rela- 
tive velocity or acceleration of the water or air with respect to 
the component that is being referenced, A is the vector of pro- 
jected areas in the x- and z-directions, V is the submerged vol- 
ume, Cd(Re) is the vector of coefficients of drag in the x-, y- 9 
and z-directions, C a is the vector of coefficients of added mass 
in the x- 9 y- 9 and z-directions, and / is the combined drag and 
added mass force vector in the x- 9 y- 9 and z-directions. The drag 
force portion of this equation is a combination of skin friction 
and residual drag. 

For the spar, skin drag is modeled only for a longitudinal 
flow along its length (parallel to the z-axis) and for the catama- 
ran, skin drag is modeled only for a longitudinal flow along their 
lengths (parallel to the x-axis). The skin drag coefficient used for 
the spar and catamaran is calculated using the Reynolds number 
calculated along the length of the spar and pontoons respectively. 
The skin drag coefficients used for the longitudinal flows, C^ n 9 
are calculated as suggested by [17], 


C^ n 


0.0075 

(l°g 10 (fle)-2) 2 ‘ 


( 2 ) 


The skin drag force varies nonlinearly along the length of the spar 
as the relative water velocity varies along the spar and catamaran 
pontoons lengths and drag force is a nonlinear function of this 
relative velocity. For this reason, numerical integration is used 
to calculate the net skin drag force along the spar and catamaran 
pontoons. 

Residual drag is calculated on the spar body, spar thrusters, 
and catamaran pontoons. For the spar, residual drag is calcu- 
lated for a flow across its body, on the bottom end (bow), and on 
its thrusters. For the catamaran, residual drag is calculated for 
a flow perpendicular to the pontoons and along the pontoons at 
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their ends. Strip theory is used to calculate the forces perpen- 
dicular to the spar and pontoon buoys while a constant relative 
water velocity is assumed for the ends of the spar and catama- 
ran pontoons. The residual drag coefficient for a flow across the 
spar (jc-direction) is assumed to have a constant value of 1 .0, as 
suggested by [18]. This same coefficient is used to calculate the 
residual drag across the pontoons (z-direction) as both the pon- 
toons and the spar are approximately cylindrical. The coefficient 
or residual drag for a flow along the spar (z-direction) at the sub- 
merged end is also assumed to have a value of 1 .0, as suggested 
by [18]. The same coefficient is used for the residual drag on 
the ends of the pontoons. Residual drag forces on the submerged 
thrusters are calculated independent from the rest of the spar as 
they are the only components of the RDSP that protrude out from 
the spar’s frame. The residual drag force on a thruster is cal- 
culated using separate but constant coefficients for flow parallel 
and perpendicular to the thrust axis. These coefficients are de- 
rived from tow tank tests, where the thrusters are towed facing 
parallel and perpendicular to the tow direction. The total drag 
coefficients are derived from fitting curves to the data, which are 
a function of velocity squared. These coefficients are used to cal- 
culate the drag forces for a flow parallel to each thruster and for 
a flow perpendicular to each thruster. 

Added mass forces for both the spar and catamaran are mod- 
eled as an increase to the virtual mass and moment of inertia of 
the RDSR This simplification can be made since waves are not 
considered for this application and therefore the local RDSP ac- 
celeration is the local relative acceleration. When operating in 
VOM, the virtual mass and virtual moment of inertia are calcu- 
lated for the RDSP in each direction using the actual mass and 
moment of inertia (including entrained water), the added mass 
and moment of inertia from the spar, and the added mass and 
moment of inertia from the catamaran . The coefficient of added 
mass for spar acceleration in the jc-direction, perpendicular to the 
spar, is chosen as 1 and the coefficient of added mass for spar 
acceleration in the z-direction, parallel to the spar, is chosen as 
1 /tt as suggested by [18]. These same coefficients are used for 
a flow perpendicular and parallel to the catamaran pontoons, re- 
spectively. 


Thruster Forces The individual thruster forces and 
thrust rate are limited to match the maximum thrust to the power 
plant suggested by [19] using Froude scaling. This limits the 
maximum thrust per engine to 47 N and the maximum thrust rate 
of change per engine to 37.2 N/s for the motor configuration cho- 
sen for the 1 / 1 0 th scale RDSP prototype with thrust values scaled 
from Thrustmaster Model TH 300N azimuthing thruster [20]. 


MODEL OF VORTEX-INDUCED DISTURBANCES ON 
THE RDSP 

Vortex-induced forces on spar platforms, such as the RDSP, 
develop when eddies are shed behind the cylindrical spar element 
as it moves relative to the surrounding fluid (Fig. 2). This occurs 



Figure 2. VORTEX-INDUCED MOTION ON A CIRCULAR CROSS- 
SECTION. 

when the relative velocity vector of the fluid has a component 
that acts parallel to the cross-section of the spar. During most 
operations, the RDSP will be dynamically positioned or moored 
in a quasi-steady current. This induces disturbance-generating 
vortices around the spar, beginning at very low relative water ve- 
locities. Vortices will also be shed at a nearly steady frequency 
and amplitude when the RDSP is traveling in VOM at a nearly 
constant speed relative to the water velocity, such as when trans- 
ferring cargo from a slowly moving vessel. The frequency of 
these eddies can be estimated, but precise knowledge of the vor- 
tex shedding frequency cannot be determined analytically. 

Eddy Shedding Frequency 

For a spar with a given diameter, the eddy shedding fre- 
quency can be approximated using the relative water velocity [5]. 
The frequency of shedding vortices from one side (see Fig. 2), 
©d = 2nfd = §?, is also the frequency of the oscillatory motion 
induced in the transverse direction. This is typically estimated 
empirically as a function of the Strouhal number with a relation- 
ship defined by [1] as 


where S is the Strouhal number, D is the diameter, and Vo is the 
relative free stream velocity. The Strouhal number itself varies 
according to a poorly defined, nonlinear function of the current 
velocity. A relationship given by [1] shows that the Strouhal 
number ranges from 0.18 to 0.32. However, it can be approx- 
imated as a constant with a value of 0.21 for cylindrical bodies, 
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as suggested by [5]. Thus, the eddy shedding frequency for the 
RDSP is modeled as a linear function of the cross flow velocity. 
For the full scale version, initial sizing estimates are in the range 
of 100 m length, 10 m diameter. Using these sizes as a reference 
and relative water velocities from 0.1 to 2.0 m/s, eddies will be 
shed at frequencies from 0.013 to 0.251 rad/s (period from 500 
to 24 s). 

Because of the Strouhal number variability and other un- 
certainties (such as unknown or unmodeled dynamics or sensor 
resolution/calibration errors), a vortex shedding frequency esti- 
mate will likely vary by as much as 15 to 20% from its actual 
value. This crude estimate of the frequency provides a strong 
motivation for using the PLL to better-estimate the disturbance 
frequency in the loop, as proposed in [7]. 

Eddy Shedding Force 

The disturbances generated by vortices are nearly periodic 
and somewhat theoretically predictable [2]. This eddy shedding 
force is calculated as [ 21 ] 

F L (f) =^z.sin(ffl d f), (4) 

where Ai is the lateral eddy shedding force amplitude, cod is the 
eddy shedding frequency that varies as a function of velocity ac- 
cording to (3), and t is the time in seconds. The amplitude of the 
force is defined as [ 1 ] 

A h = l -CaDpV%, (5) 

where L is the length of the section of interest and Cl is the co- 
efficient of lift estimated as 0.28 [ 21 ]. 

CONTROL DESIGN 

Typical control algorithms wait for error accumulate from a 
disturbance before attempting to compensate for it, but the distur- 
bance rejection technique utilized here has has the unique ability 
to anticipate and reject disturbances before they induce motions 
on a system. This pre-existing adaptive control and disturbance 
rejection algorithm [10] is augmented with a Phase Locked Loop 
for frequency synthesis. The result of this inclusion is improved 
performance and robustness with respect to variations in distur- 
bance frequency. 

Direct Adaptive Output Regulation And Disturbance 
Rejection 

Model Reference Adaptive Control (MRAC) with Direct 
Adaptive Disturbance Rejection (DADR) has been shown to mit- 
igate persistent disturbances of known frequency but unknown 


amplitude and phase [10,22-26] by taking advantage of periodic 
or persistent disturbance behavior and adaptively estimating the 
gains. When applied to large space structures, ocean platforms, 
wind turbines, and evolving systems [11, 14,27-30], this con- 
troller has produced good results. 

The theoretical development of the controller assumes a Lin- 
ear Time-Invariant (LTI) plant with the addition of a persistent 
disturbance, given by 

Xp = ApXp -|- BpUp -|- I> d 

y p = CpX p ; x p (0)=xjj, 

where x p is a vector of plant states, y p is the output, w d is a per- 
sistent disturbance, and u p is the control command. 

When applying the adaptive control algorithm from [10] to 
a regulation problem, the goal reduces to defining a control law 
and adaptively adjusting the controller gains such that the plant 
output, y p , given by ( 6 ) will go to zero. This is done in conjunc- 
tion with adaptive rejection of persistent disturbances. For regu- 
lation, the reference system input and state become zero. Thus, 
the adaptive control law from [ 10 ] reduces to 

Up = Gyp + H$<u (7) 

where <|) d is a vector that forms a basis for the disturbance. The 
corresponding adaptive parameter update law given in [ 10 ] re- 
duces to 

G = A\y P yl~aG, 

T (o) 

H= A4y p (|> d - a H, 

where Ai, A 4 , and a are design parameters. A block diagram of 
this control setup is shown in Fig. 3. 



Figure 3. BLOCK DIAGRAM FOR ADAPTIVE OUTPUT REGULATION 
WITH DIRECT ADAPTIVE DISTURBANCE REJECTION. 
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Adaptation For Poorly Known Frequencies 

When the frequency of a disturbance is not precisely known, 
a complete basis can no longer be predefined for this disturbance. 
In this case, the performance of the DADR scheme is degraded 
[11, 14], Here, a Phase Locked Loop (PLL) is used to better- 
estimate the disturbance frequency. This allows the functions in 
(|)d within the DADR algorithm to more closely approximate the 
true basis. The focus of this modification is for the rejection of 
a single vortex-induced disturbance, but the methodology is also 
applicable for plants with multiple disturbances. 

Phase Locked Loop A Phase Locked Loop (PLL) is a 
circuit that synchronizes the frequency and phase of an output 
signal to that of a given input signal [31, 32]. The PLL that is 
used here is shown in Fig. 4. Boxes highlight the three main 
components that are common to all PLLs: synchronized oscil- 
lator, phase detector, and loop filter. In general, the role of the 


PHASE DETECTOR 



Figure 4. PHASE LOCKED LOOP DETAILS. 


three main components may be described as follows. The phase 
detector component of the PLL compares the input and the Syn- 
chronized Oscillator signals to determine an error signal that is 
proportional to the phase difference between the two. This er- 
ror feeds through the Loop Filter (LF) and drives the frequency 
in the opposite direction. By continually updating the frequency 
in the Synchronized Oscillator in this way, the signal is better- 
synchronized with the input signal in terms of phase and fre- 
quency. The PLL is designed to reduce the phase error until the 
output of the LF becomes stationary or “locked”, thus inspiring 
the name “Phase Locked Loop”. When operating in this locked 
state, the frequency exiting the LF will match that of the input 
signal. 

The particular formulation shown in Fig. 4 is a simplified 
version of [33], and can be found in [34—36] for the case where 
G(s) = G. Through linear analysis of the PLL, G(s) from [34,35] 
is re-defined as G(s) = to provide an additional zero at 


the origin of the closed loop transfer function. Its purpose is to 
cancel with the extra pole generated by a disturbance input with 
a step change in the frequency rather than in the phase. This 
modification is made to improve the tracking performance. 


Modified Control Framework The control framework 
set forth in [10] is modified to include an on-line frequency esti- 
mate that is provided by a PLL. A block diagram of the adaptive 
controller and DADR with the inclusion of a PLL is shown in 
Fig. 5. The RDSP’s output goes through a bandpass pre-filter 
that extracts the component of the signal within the frequency 
range of interest. This filtered output is used by the PLL to refine 
the estimate of the disturbance frequency. The PLL frequency 
estimate is lowpass filtered to eliminate any high frequency os- 
cillations. Finally, this frequency replaces the static frequency 
guess in the DADR basis functions. 



Figure 5. CLOSED LOOP SYSTEM WITH M RAC/DAD R AND A PLL. 


The bandpass pre-filter allows all of the signals outside the 
range of interest to be removed. Therefore, the input to the PLL 
primarily contains the component of the plant’s output that is a 
direct result of the disturbance. If multiple disturbances are ex- 
pected, the framework may be adapted by including additional 
sets of PLLs with corresponding filters. For example, n distur- 
bance frequency estimates would require n sets of bandpass fil- 
ters, PLLs, and low pass filters. Each bandpass pre-filter will be 
adjusted according to the expected frequency of the correspond- 
ing disturbance. 

The filtered PLL frequency, C 0 p£ L is used in the DADR 
scheme rather than the unfiltered frequency, copll, to reduce the 
oscillations in the DADR basis functions. This also reduces the 
effect of the disturbance rejection portion of the controller at the 
beginning of the simulation, giving the PLL time to lock onto the 
desired frequency without having a large amount of interaction 
with the DADR. 
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RESULTS 

For each 6 -DOF simulation, the six forward facing pro- 
pellers provide a constant thrust of 300 N (50 N each) in the 
x-direction, inducing a surge velocity of approximately 0.22 m/s. 
An initial surge velocity of 0.22 m/s is given to the RDSP but 
there is still a transient period before the disturbance frequency 
converges to a constant value (Fig. 7, Fig. 11). A basic PID 
controller is used to minimize the RDSP’s yaw motion while the 
adaptive controller is configured to reduce the transverse oscilla- 
tory motions. As with the 1 DOF model [7, 14], the control law 
given by (7) with parameter update laws given by ( 8 ) is used in 
conjunction with a PLL. This controller implementation is shown 
in Fig. 5. 

The gains in the update law ( 8 ) are chosen to be Ai = 
200 , 000 / 2 , A 4 = 20 , 000 / 2 , and a = 1 ■ 10~ 4 . For the PLL, the 
gains are chosen to yield negative (stable) roots for s 3 + as 2 + 
aGs + aGKi = 0 by selecting G = 2.5 • 10“ 4 , K t = 2.844 ■ 10~ 6 , 
and a = 0.3. 

The frequency guess for vortex disturbances was taken as 
c os = 0.25 rad/s. This is approximately 78% of the steady-state 
frequency of vortex-induced disturbances (cod ~ 0.32). Thus, 
|®d — o>s| is selected such that it exceeds the expected frequency 
estimation error for this application. 

For this simulation, two sensor/actuator configurations are 
considered, as shown in Fig. 6 . The first is a Single-Input-Single 


a) SISO a) MIMO 




Figure 6. SENSOR AND ACTUATOR LOCATIONS. 


Output (SISO) scenario, where the sensor is located at a point on 
the catamaran (near the top of the spar) to detect motion in the 
y-direction. The control input is taken to be a point force acting 
at the location of the central thruster, or equivalently equal thrust 
from all lateral thrusters. For the second simulated configuration, 
the input and output are collocated at the center thrusters. In this 


case, the output contains both the sway and roll velocities and 
the input is the force and moment in these directions. These out- 
puts can be estimated by measuring the motions at the catamaran, 
where the sensor system is currently located on the RDSP proto- 
type, and transforming these measurements to motion estimates 
at the center set of thrusters. The force input can be applied to 
the physical system by dividing the total desired force equally 
among each of the lateral thrusters, as with the SISO simulation. 
The moment portion of the input can be achieved by superim- 
posing equal and opposite thrusts from the top and bottom sets 
of lateral thrusters. 

The effectiveness of the controller when used in conjunction 
with an on-line frequency estimate was shown for a simple 1 - 
DOF simulation in [7, 14] and will now be verified for a 6 -DOF 
RDSP model with increased dynamical complexity. 

Single-Input-Single-Output Scenario 

For the SISO scenario, the output is taken to be the velocity 
of the 1 / 10 ** scaled prototype in the y-direction at the top of 
the spar. The input is the force in the y-direction, applied at the 
middle thrusters (or equally across each of the lateral thrusters). 
This simple sensor/actuator configuration may be able to realize 
the goal of mitigating the vortex-induced transverse oscillatory 
motion at the top of the spar, but the control system is not capable 
of intentionally eliminating RDSP pendulations about this point 
since they are not visible in the output (see Fig. 6 ). 

The on-line frequency estimate is shown in Fig. 7. In Fig. 



Figure 7. DISTURBANCE FREQUENCY ESTIMATE FOR THE SISO 
SIMULATION. 


8 , the disturbance is shown alongside the DADR portion of the 
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controller at the output (sensor) location. The standard devia- 
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Figure 8. DISTURBANCE REJECTION EFFECTIVENESS FOR THE 
SISO SIMULATION: FULL RUN TIME (LEFT) AND LAST 100 s (RIGHT). 



Figure 9. SWAY VELOCITY COMPARISON FOR THE SISO SIMULA- 
TION. 


tion of the acceleration at this point due to the vortex-induced 
disturbance, given by C p r p «d is 1.623 ■ 10 -3 m/s 2 . When the 
DADR with a PLL is included, the STD of (CpT p u^ + C p £ p //c|>d) 
becomes 1 .056 • 1 0 -3 m/s 2 , and thus reduces the disturbance ap- 
pearing in the output by 35%. 

The portion of the error not rejected a priori by the DADR 
scheme enters into the dynamical system (RDSP) and must then 
be mitigated by the output feedback term in the adaptive control 
law. The open loop STD of the velocity resulting from the vortex 
shedding for the last 1000 s of the simulation is 3.046 ■ 10 -3 m/s. 
If the loop is closed with the MRAC (MRAC/DADR), this is re- 
duced to 2.626 ■ 10 -3 m/s (1 .300 ■ 10 -4 m/s). The inclusion of the 
PLL with the MRAC/DADR algorithm yields a transverse veloc- 
ity error with a STD of 6.4 ■ 10 -5 . Thus, the controller nearly 
eliminates the velocity error, y p , at the top of the spar (as shown 
in Fig. 9). After removing the bias, the open and closed loop 
displacements are shown in Fig. 10 for the latter 1000 s of the 
simulation. As with the velocity, the oscillations in position are 
almost completely eliminated. 



Figure 10. TRANSVERSE DISPLACEMENT COMPARISON FOR THE 
SISO SIMULATION. 


Multi-Input-Multi-Output Scenario 

The output for the MIMO example is a vector containing the 
sway velocity in m/s and 10 times the roll velocity in rad/s. The 
sway velocity output signal is selected as the input to the PLL, 
which provides an estimate for the frequency of both the force 
and torque disturbances. Because the disturbance frequency it- 
self oscillates about an equilibrium, the PLL gains were selected 
to provide a somewhat sluggish response. Thus, it takes 600 s 
for the frequency estimate to initially approach the actual distur- 
bance frequency (as with the SISO case). This is shown in Fig. 
11 . 


Fig. 12 shows the disturbance input («d) plotted against the 
DADR component of the control command (Hi |>d) with and with- 
out the inclusion of a PLL. This is shown over only a portion of 
the simulation time in order to provide a closer look at these sig- 
nals. The standard deviation of the disturbance force and mo- 
ment over the last 1000 s of the simulation are 49.82 N and 
36.05 N-m, respectively. The disturbance rejection component 
of the control command, Hi j>d, does a much better job at oppos- 
ing the disturbance signal, «d> than in the SISO scenario. This 
may also be observed statistically by again considering the unre- 
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Figure 1 1 . DISTURBANCE FREQUENCY ESTIMATE FOR THE MIMO 
SIMULATION. 
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Figure 12. EFFECTIVENESS OF DISTURBANCE REJECTION FOR 
THE MIMO SIMULATION. 

jected portion of the disturbance. Without the inclusion of a PLL, 
STD(wd +ZJ(|>d) is 2.55 N (95% reduction) for the sway distur- 
bance component and 10.54 N-m (71% reduction) for the roll 
disturbance component. When these STD values are calculated 
again with the inclusion of a PLL, the result is 1 .47 N (97% re- 
duction) for the sway disturbance and 8.55 N-m (76% reduction) 
for the roll disturbance. Thus, approximately three times as much 
of the disturbance is rejected a priori for this sensor/actuator sce- 
nario (and, in fact, almost all of it). 

As could be expected based on the effectiveness of the 
DADR component alone, the controller results in a significant 


reduction in output error (sway and roll velocities) in all cases. 
For the open loop results, the sway velocity has a standard de- 
viation of 7.1 ■ 10~ 3 m/s and the standard deviation of the roll 
velocity is 3.88 ■ 10 -2 deg/s at the center thruster over the last 
1000 s. When the MRAC/DADR with PLL is applied, the STD 
is reduced to 1.3 1 ■ 1 0 -4 m/s for the sway velocity and 8.36-1 0 -3 
deg/s for the roll velocity. This corresponds to a reduction in er- 
ror of 98.2% for sway velocity, v - more than a two-fold increase 
in overall output error reduction when compared to the SISO re- 
sults. The roll velocity, p 9 was similarly reduced by 95.2% (this 
DOF was not controlled in the SISO case). The corresponding 
displacements in the sway direction are shown in Fig. 13, and 
the roll angles are shown in Fig. 14. 



Figure 13. TRANSVERSE DISPLACEMENT COMPARISON FOR THE 
MIMO SIMULATION. 


CONCLUSIONS 

An adaptive control and disturbance rejection algorithm 
with the inclusion of disturbance frequency synthesis was shown 
to mitigate vortex-induced motions for the RDSP operating as a 
spar platform in VOM. Simulation results were shown for two 6- 
DOF RDSP actuator/sensor scenarios: (i) a single control force 
located at center thruster and a single output at the top of the spar, 
and (ii) control actuation using both force and moment generated 
at the center thruster with two collocated outputs. The location 
of the actuators and sensors in the latter (MIMO) scenario re- 
duced the transverse oscillitory motion to half of that found for 
the former (SISO) scenario. In both cases, the inclusion of a PLL 
allowed the MRAC/DADR to more effectively mitigate the per- 
sistent, vortex-induced disturbance at a frequency that was not 
well known. 
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Figure 14. ROLL ANGLE COMPARISON FOR THE MIMO SIMULA- 
TION. 
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